Non-conditioned generation of Schrodinger cat states in a cavity 
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We investigate the dynamics of a two-level atom in a cavity filled with a nonlinear medium. 
We show that the atom-field detuning S and the nonlinear parameter x rnay be combined to 
yield a periodic dynamics, allowing the generation of almost exact superpositions of coherent states 
(Schrodinger cats). By analysing the atomic inversion and the field purity, we verify that any initial 
atom-field state is recovered at each revival time, and that a coherent field interacting with an 
excited atom evolves to a superposition of coherent states at each collapse time. We show that a 
mixed field state (statistical mixture of two coherent states) evolves towards an almost pure pure 
field state as well (Schrodinger cat). We discuss the validity of those results by using the field fidelity 
and the Wigner function. 
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INTRODUCTION 



In quantum optics, the Jaynes-Cummings model (JCM) de- 
scribes the interaction between a two-level atom and a single 
quantized mode of the radiation field in a lossless cavity and 
within the rotating wave approximation (RWA). The JCM is 
probably the simplest fundamental model of field-matter inter- 
action with an exactly integrable Hamiltonian. Just over forty 
years since its introduction P|, the model has originated sev- 
eral studies in various contexts and with different purposes, and 
has become the basis for several generalizations and other models 
Q]. More recently, important experimental achievements in cavity 
quantum electrodynamics (QED) and trapped ions have stimu- 
lated both theoretical and experimental research in that area Q. 
An interesting related subject is the quest for generation meth- 
ods of macroscopically distinguishable superpositions of quantum 
states, or Schrodinger cat states H. Several schemes using co- 
herent states have been proposed |M |y, Q, |M l3> 8'iid a few ex- 
perimental realizations have been already accomplished in cavity 
QED as well as in trapped ions systems llOl . [ill 11^ . In cavity 
QED models, states close to those superpositions arise at specific 
times, for the cavity field initially in a coherent state p^ or even 
in a statistical mixture of two coherent states \iM . Propositions 
such as the Yurke-Stoler generation scheme p\ , and those based 
on quantum non-demolition processes (QND), depend on very 
large values of Kerr nonlinearities, which is probably the main 
obstacle for their implementation ilOi] . However, in the last few 
years, the observation of large Kerr nonlinearities with low inten- 
sity light [la . Ha. 1171 as well as propositions involving small Kerr 
nonlinearities [lal have renewed the interest on those schemes. 
Furthermore, schemes for generation involving cavity QED with 
a nonlinear medium, based on atomic conditional measurements 
have also been proposed |19|. 

In this paper, we present a method that does not depend on 
conditional measurements. We have found that, the JCM with 
a nonlinear Kerr-like medium, under suitable combinations of the 
atom-field detuning S and the nonlinear parameter x and for an 
initial field prepared either in a coherent state or in a statistical 
mixture of two coherent states, makes possible a Schrodinger cat 
state generation with higher fidelity than the JCM without a 
nonlinear medium. We would like to remark that in the ordinary 
JCM, the initial field in a coherent state evolves to a state close to 



a Schrodinger cat state at a specific time, as reported in jTj], and 
if we start with a statistical mixture of two coherent states, the 
field basically remains in a mixed state p^l2(l]. The possibility of 
generating superpositions of coherent states in the JCM with a 
nonlinear Kerr-like medium with an atom-field detuning has not 
been yet addressed in the literature. 

This paper is organized as follows: in section II we introduce 
the model and obtain the evolution operator in the RWA approx- 
imation. In section III we present the numerical results of some 
fundamental quantities and show how to obtain the condition for 
a periodic dynamics. In section IV, we discuss the main results 
and present our conclusions. 



II. MODEL 

In this section we describe the interaction of a two-level atom 
with a high-Q single-mode cavity filled with a nonlinear Kerr- 
like medium, which can be modelled as an anharmonic oscillator 
|21l . The cavity field is coupled with both the two- level atom 
and the nonlinear medium. If the response time of the nonlinear 
medium is sufficiently small we can adiabatically eliminate the 
photon-photon coupling, i.e. considering the field and nonlinear 
medium frequencies far from each other 1221 . Then, the total 
Hamiltonian of the system, with the adiabatic, RWA and dipole 
approximations, can be written as |24l | 



H = fiLUQa a - 



^huj^gCrz + fix a a + hQ,(a(T_ + acr^), (1) 



where loq (u)eg) is the cavity field (atomic transition) frequency, 
a^ (a) is the creation (annihilation) operator of the cavity mode 
obeying [a,a'^] = 1, a^ = \e)(e\ - \g)(g\, fT+ = \e)(g\ and 
<T_ = \g)(e\ are the standard Pauli matrices operators, where |e) 
(\g)) refer to the excited (ground) atom state, Q, is the atom-field 
coupling constant and x is the nonlinear parameter, propor- 
tional to the dispersive part of the third-order nonlinear suscep- 
tibility 113. 

Following the approach of Stenholrn ze|, we delineate the main 
steps to obtain the exact (under the RWA) time evolution opera- 
tor for this model. After some algebra, we can rewrite equation 1^ 
as 



H — Hq + if in 



(2) 
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where 



hLOQ(a^a+ IfT^), 



-Hint = -ff 1 + -ff 2 , 



(3a) 
(3b) 
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with 



Hence, the evolved atom-field state is given by 



fl"i = ^hx,^-" + hx^-"[{a^ ay + a^ acr,] - hx^-"{a^a + ^a,), (4a) 

H2 = fi[| - x'^Ha^o - h)]<^. - h^^X^'^^ + fin(a^o-- + aa+), (4b) 

where 5 = tJeg — ajp is the atom-field detuning. 

We have verified that equation l|^) and equation l^j) com- 
mute, and, therefore we may write Hj — Uo{t)H^^JjQ{t) — i/int, 
which is just the Hamiltonian in the interaction picture. Hence, 
the respective time evolution operator is given by 

Ui{t) = C/i(f)C/2(t) = exp {-jH^t) exp {-^H^t) , (5) 

where the exponentials have been decoupled. After some manip- 
ulation we obtain the following form 



with 



TT l-t\ — [ -^n + l ^ \ ( ^n + 1 -B„_|_i 



-Bji+i — e 



A ^lo +^ sin(in„+it) 
A„+i =cos(2n„+it) -i7„+i j4 , (7b) 



f^n + l 



Bn+-i ~ 2ii}a 



■|. sin(i»„_|_it) 

f^n + l ' 



where 



f^n+i = \/7^+i+4n2(n + l), 



(8) 



and7„+i = 5 - 2x*^'n. 

The result in equation Q has been already obtained, in another 
context Ull. 



III. ATOM-FIELD DYNAMICS 



In what follows we are going to assume an uncorrelated initial 
atom-field state, i.e. 



P = Pa'XJPf, 



(9) 



where p^ is the initial atom density operator, initially an ex- 
cited state* Pa = |e)(e| and pi is the initial field density op- 
erator, initially either a coherent state Pf'^ — \a}{a\ or an 
equally weighted statistical mixture of two coherent states pf^ = 
|(|a)(a| + | — a)(— a|). In all cases a = |Q:|e"^, and we will fix 
\a\ = y/n — 5. The general form of the initial field state in the 
Fock state basis is 



Pc = '^Pn.,n\n){rn\, 



(10) 



where Pn,m = {'n\p[\m} are the initial field matrix elements. For 
the coherent state 



Pn.n 



-\a[ 



\/n\m\ 
and for the statistical mixture of two coherent states 



(11) 



p{t) = Uj{t)pUl{t) = 



\l^\ _ ( Pee{i) Peg{t) 



(13) 



Pg^it) Pggit) 

whose elements in the atomic basis are 

Pee{t) = ^ /5n,m-E„+i-E;^+i^„+iA^+iin)(mj, (14a) 

Pegit) = -^ Pn,m.En+iE^+iAi+iB^+i\n}{rn + 1|, (14b) 

n,77i 
Pgeit) = - y~] Pn,mE„ + iE:^ + iB„+lA*m+l\n + 1) (m] , (14c) 

Pggi^-) = "^ Pn,mE„+iE'^^iB„+iB^_f_i\n + l){m+ 1\, (14d) 



(6) with 



(7a) 



E„+i — e 



(15a) 



,1 sin (ir^n+if) 
A„+i =cos(^n„+if) -J7„+i ^ , (15b) 



B„+i = 2in^^rrT 



sin(it7„+it) 



"„+i 



(7c) being 



"n+i = ^7^+1 + 402 (n + 1), 



(15c) 



(16) 



the generalized Rabi frequency, with j„+i = S — 2x n. 

Once that p{t) belongs to the trace class operators acting in 
the space corresponding to the direct product in equation @ we 
can trace over the field variables in equation 1131 to obtain the 
reduced atomic density operator 

V -^eg ^gg J 

where A^^ — '^{n\pij{t)\n) . Using equation CU, we have 

n 

Pa(t) = EPn,nlAi + lHe)(el + E Pn,n 1S„+1 ^ Ig) (g] 
n n 

-E(pn + l,ne-'^*"('" + l)*A„ + 2B: + i|e)(g| +C.C). (18) 
n 

Analogously, by tracing over the atomic variables, we obtain the 
(reduced) field density operator 



Pfit) = Trjp(i)] = ^p„^„(i)|n)(m|. 



(19) 



where 



Pn,mii) = {'rAPf(t)\-m) = En^^El^^^^Pn^^A^+iA*^^^ 

I Pn— l,m — 1^ ^-'n^-^mji 

are the evolved field matrix elements. 



A. Atomic Inversion 



(20) 



Pn 



sm 1 — |a| 



m 2 



Vnlmi 



[l-f(-l)"+'"]- 



^ ' A quantity usually measured in experimental cavity QED is 

the atomic population inversion |2a . 12911 , defined as the difference 
between the probabilities of finding the atom in the excited state 
and in the ground state. Here the atomic inversion is given by 



* The extension to a more general initial atomic state, like p^ 

e 

E Pi,j\i){j\ where p^j = {i\pjj), may be easily done. 



W(t) = Tr,[o-,p,(t)] = Y, Pni\Ar,+ if - 1B„+ 



11^), 



(21) 
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FIG. 1: Atomic inversion as a function of t/t^ when the field is 
initially in pf = |a){a| or pf^ = |(|a){a| + |-a)(-Q:|) (q = 5, 
= 0) and the atom is initially in p^ = |e)(e| with 5 = and 

X^^' = 0. 



FIG. 2: Atomic inversion as a function of t/t^ when the field is 
initially in pf = \a){a\ or pf^ = i(|a){a| + |-a)(-a|) (a = 5, 
= 0) and the atom is initially in p^ = |e)(e| with 5 = S^ = 4.8f2 
and x*^' = O.m. 



where P„ = p„ „ = {n|pf|n) is the initial field photon number 
distribution. 

It is well known that the atomic inversion is very sensitive to 
the initial field photon number distribution P„. For the field 



in the Fock state p\ = |n){n| we have P„ 



resulting a 



sinusoidal behaviour for VV(i). If the initial field is the thermal 
state pf" = J2 -Prfl^) {n\ we have P,f = n"/{n+l)"+^ , and a more 

n 

irregular behaviour occurs 29]. For the field states considered in 

this paper we have P^'' = P™ = e^'"' JQJ^"/n! so that the atomic 
response, at least regarding the atomic inversion, is the same in 
either case. The atomic inversion reveals non-classical features: 
the Rabi frequency oscillations present collapses and revivals [30|. 
We have parametrized the variable time as t/t^ to allow a better 
comparison among the different plots in terms of the revival time 
t^. In figure we plot the atomic inversion as a function of t/t^, 



having S = x 



(3) 



0, and we observe the pattern of oscillations 



characteristic of the ordinary JCM atomic dynamics. 



B. Linear Rabi Frequency 

We would like to find under which circumstances we may have 
a periodic dynamics. One way of doing that is to treat the Rabi 
frequency as a continuous quantity, so that we may expand equa- 
tion 1161 1 around the initial mean photon immber n 



^r, 



E 



1 d''n„ 



k\ 



in- 



(22) 



periodic behaviou r l3ll. This is the case, e.g. for the intensity- 
dependent JCM p^. l32l . We show that it also may be the case 
for the JCM with a Kerr-like medium: from the second order 
derivative of the Rabi frequency. 



f^„ 



(3)2q2 



A^ 



"n + l 



(24) 



we have H^^i = if A^^i = x'^^'^n+i or equivalently 

-2x'-^\ (for all n). (25) 



S^-^^-O.M^ 



2x(3) 



It is clear from equation (I22II that all higher-order derivatives 
up to the first-order vanish when 5 — 5^. This condition deter- 
mines the periodic behaviour in the dynamics of the model. As a 
first illustration of that, we plot in figure |5| the atomic inversion 
for the same conditions of figure but satisfying the relation for 
5^ in equation 1251 above which assures the periodic behaviour. 
Furthermore, if we insert equation 12511 in equation 12311 . we ob- 
tain 









and equation llbll becomes 



n. 



:(5, + 2x''>(7i + 2). 



(26) 



(27) 



The first term above governs the rapid oscillations in the Rabi 
frequency while the remaining terms generates the envelopes (re- 
vivals, super-revivals and so forth) . It is well-know that two suc- 
cessive terms (in the discrete spectrum) of Rabi frequency, i.e. 
n^_i_i and Qfi have a 2-k phase difference, so that the revival time 
is given by 



27rn: 



^71 + 1 



^ri+1 



(23) 



where A^+i = Q, — x 7ri+i- If only the first two terms in 
equation 12211 are nonzero, the Rabi frequency exhibits a perfectly 



We remark that similar results were obtained, in another con- 
text, in the two-photon JCM with Kerr-like medium [2511 and in 
|3al . where the authors obtained the linearized Rabi frequency, 
although they just discussed the behaviour of the atomic inver- 
sion and used a strong-field approximation {n^ 2> n) to obtain 
the evolution operator. 

We would like now to comment about the physical relevance of 
the values of S taken in this paper, i.e. if they are consistent with 
the RWA. From experimental realizations in microwave cavity 



QED [THHIll, we have that Q ~ lO^'Hz, 
loq ~ lO^^Hz. Here we are considering S -^ 
consistent with the RWA once that 6 ~ 10"'* 



igg ~ 10 Hz, and 
lO^fiHz which is 
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FIG. 3: Field Purity as a function of t/tr when the field is initially 
in (I) pT = \a){a\ or (II) pr = i(|a)(a| + \-a){-a\) {a = 5, 
= 0) and the atom is initially in p^ = |e)(e| with 5 = and 
X^^' = 0. 



C. Field Purity 

A very useful operational measure of the field state purity is 
given by the linear entropy 



Cf(t) = 1 - Trj [pf (i)] = 1 - ^ |p„,™(t)| 



(28) 



In figure |21 we plot the linear entropy for the resonant case 
and in the absence of the Kerr-like medium, i.e. with x = 0. 
It is well-known [la| that at half of the revival time (collapse 
region), the initial coherent field evolves towards a field close to 
a pure {Schrodinger cat) state, figure |S}I, whereas for an initial 
statistical mixture of coherent states, the field is always far from 
a pure state 20], as shown in figure EJII. 

The situation is very different if we consider the condition that 
gives a periodic dynamics. In figure 2] we plot the field purity 
for the same conditions as considered in the calculation of the 
atomic inversion. The initial coherent state evolves to an almost 
pure state (very close to a superposition of two coherent states, as 
we will discuss in what follows) at each collapse time and returns 
to the initial state at each revival time, figure|lj-l. Remarkably, an 
initial statistical mixture of two coherent states also evolves to an 
almost pure state (approximately a superposition of two coherent 
states) at each collapse time and returns very close to the initial 
state at each revival time. As seen in the atomic inversion plot 
in figure l^l the initial atomic state is basically recovered at each 
revival time. 

It is well known [3, 12^ that a nonlinear Kerr- like medium may 
convert a field in a coherent state to another pure state, namely 
the Yurke-Stoler Schrodinger cat state. In the model present here, 
however, a field in an incoherent superposition (mixed mixture) 
can, under suitables conditions, evolve to a state very close to a 
pure state (a coherent superposition of two coherent states). In 
this case of course the other subsystem (atom) is left in a mixed 
state. 
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FIG. 4: Linear entropy as a function of t/t^ when the field is 
initially in (I) pf = |a)(Q| or (II) pr = |(|«){q| + |-a){-a|) 
{a = 5, (j) = 0) and the atom is initially in p^ = |e){e| with 
S = 5^ = 4.80 and x*^' = 0.1!^. 



quasi-probability distribution which is the Fourier transform of 
the anti-normally ordered quantum characteristic function |3R| . 
l3q| . For the field calculated here, the Q-function is given by 



0(/3,t) = -(/3|pf(i)|/3) 
n 



-W 



E- 



M 



\Jn\m\ 



(29) 



where ]/3) is a coherent state with /3 = 5R(/3) -I- ^^(/S). 

In discussions found in the literature, the Q-function is plotted 
only at some specific times |2(ll25|l. as shown in figure|3 In order 
to provide more complete information about the field evolution, 
specially at the collapse time, when the pure state generation oc- 
curs, we present the full time evolution of the field Q-function as 
a form of an animation. The animation, external to this paper, 
can be downloaded^ from here Interesting results arise when the 
dynamics is periodic: there is a relation between fractions of the 
revival time and the number of peaks in which the initial coher- 
ent field splits (or recombines), as shown in figure|S] For instance. 



D. Q and Wigner Functions 

In this subsection we consider the field dynamics from the point 
of view of the Q and the Wigner functions. The Q-function is a 



t It is necessary to download the file from http:// 
j ournalsonline . tandf . co . uk/openurl . asp?genre=article&id=doi : 
10.1080/09500340500058116 and follow the instructions therein. 




FIG. 5: Q-function at (I) t = 0, (II) t = \t^, and (III) t = \t^ 
when the field is initially in {a,c) pf — \a){a\ (a = 5, = 0) or 
{b,d) p™ = i(|a){a| + \-a){-a\) (a = 5, = 0) and the atom 
is initially in p^ = |e)(e| with {a,b) <5 = and x ~ and (c, d) 
5^5^= 4.8r2 and x*^' = O.m. 



at half-revival (collapse) time, the field becomes basically a pure 
{Schrodinger cat) state, represented by two peaks and an interfer- 
ence (oscillating) structure in phase-space. A similar behaviour 
occurs if the initial field is a statistical mixture of two coherent 
states; a pure Schrodinger cat state generation is almost perfect 
at half-revival time, although the initial state is a mixed state. 

For a better visualization of the field state generated at the col- 
lapse time we consider the Wigner function, a quasi-probability 
distribution given by the Fourier transform of the symmetrically 
ordered characteristic function |36| . Alternatively it may be writ- 
ten as [S^l 

W{l3,t) = -y^i-ir{n,l3\pi{t)\n,f3) 

n 

= -J2i-irPr.,Ut){MD{2f3)\n), (30) 

IT ^ ^ 



where Z)(2/3) — e^"^" ^ °' is the Glauber displacement operator 
and 



{m\D{2f3)\n) 



-2\l3[ 



-2|/3|- 




m~n T (m — n) 



im'- 



(m > n), 



(31) 



n-~m J- (n — m) 



im) (m ^ n). 



where L^'"^' (4|/3p) are the associated Laguerre polynomials |3a| . 

In figure |H| we plot the Wigner function at the collapse time 
for the ordinary JCM, when the field is initially in the coherent 
state and the atom is initially excited. That corresponds to a 
state close to a Schrodinger cat state, in agreement with our field 
purity analysis, figure |HJI. 

When the condition for periodicity is fulfilled, the Wigner func- 
tion at the corresponding time gives the state depicted in figure|7] 
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FIG. 6: Wigner function at f = ii^ when the field is initially in 
pp — \a){a\ {a — 5, (j) ~ 0) and the atom is initially in p^ = |e)(e| 
with 5 = and x^^^ = 0. 
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FIG. 7: Wigner function at f = ^i^ when the field is initially in 
pf — \ci){a\ {a — 5, (f) ~ 0) and the atom is initially in p^ = |e)(e| 
with 5 = 5^ = 4.8n and x^^' = O.m. 



We clearly see that the field is much closer to a Schrodinger cat 
state in this case. As we shall see, when the field is initially in 
a coherent state, all values of S^ allow an almost Schrodinger cat 
state generation at each collapse time, with each cat state hav- 
ing a specific relative phase. For complementarity, we discuss in 
appendix A the situation of large detunin g (d ispersive approxi- 
mation) based on an effective hamiltonian |39l . 

We would like to point out that the action of the Kerr-like 
medium itself is not enough to 'purify' the field. The Kerr-like 
medium merely converts the statistical mixture of two coherent 
states into another statistical mixture, as we show in appendix B. 
From the phase space point of view, the initial statistical mixture, 
figure |31-(b), splits in two deformed peaks, one for each coherent 
state of the incoherent superposition, in such a way that there is 
no perfect phase recombination at the collapse time. In the case 
of periodic dynamics, the peaks are uniform, as shown in figure|S} 
(d), and become virtually indistinguishable at the time they cross 
each other, which means that they correspond to an almost pure 
state. 



E. Mean Photon Number 

Now we present the procedure adopted to determine the rel- 
ative phases of the superposition attained at each collapse time, 
depending on the values of S and x which satisfy equation 12511 . 
The first thing to note is that the mean photon number at the 



x(3) 



n{^Q 







25.500 



0.5!r2 


25.074 



0.450 
25.110 



o.sn 
i.oen 

25.179 



0.2!r2 

2.m 

25.316 



o.ioen 

4.474160 
25.493 



O.IO 

4.80 

25.495 



0.05O 

9.90 

25.324 



O.OIO 
49.980 
25.020 



0.005O 
99.990 
25.005 



TABLE I: Mean photon number a,t t — if,, when the field is initially in p{^ — \a){ 
and the atom is initially in Pa = |e){el for difi^erent values of 5 and x • 



a\ or pr = |(|a)(a| + |-q)(-q|) (a = 5, <?!> = 0) 



r 



collapse time may not be the same as the initial one. 

To verify this we consider the mean photon number given by 



i(f) = Trf[pf(f)n]=^P„(f)r 



(32) 



where P„(f), for the model presented here, is given by 

P„{t) = /5„,„(i) = {n\pf{t)\n) = P„\A„+if + P„_i|B„|'. (33) 

As shown in table|lj for each combination of S and x we have 
a specific mean photon number at the collapse time. 



F. Field Fidelity 

To obtain the values of the relative phase tS, of the superposition 
attained at the collapse time, we calculate the field fidelity, defined 
as 



■^f(i) = Trf [pfPf(f)] = ^ pm,nPn,m{t), 



(34) 



so that the evolved field state equals the initial field state if and 
only if J-i{t) = 1. For that, we compare the field state obtained at 



J 



t — ^t^, when the condition for the periodic dynamics is satisfied, 
to the following {Schrodinger cat) state 



|5;i9) =C2(|5)+e'''h5)), 



(35) 



Pt = 



„cat 

Pf 



= |5;i9){5;i?|, where C = §(1 + e"^'"' cost?)"^ 
is the normalization constant. From the mean photon number 
values at the collapse time we have |5p = n{^tr) with a = lale'*^, 
where the value — 7r/2 was obtained from the previous analysis 
of the Q-function (see figure 1^1. We then vary the values of the 
phase i9 at the collapse time until we obtain J-{{t) « 1. The results 
are presented in table|n]for the initial coherent field state and the 
atom excited for different values of 6 and x • Apart from the first 
combination (ordinary JCM), all the others satisfy equation I25t 
and we have the generation of a state very close to a Schrodinger 
cat state at each collapse time with a specific relative phase ■&. 
We have also found that, for an initial coherent field state with 
(j> — IT, i.e. p'l^ — |— a)(— Of], the relative phase of the generated 
superposition is given by — 19, instead. We have payed special 
attention to the combinations (i) S = 4.80 and x = O.IO, which 
generates a state very close to an even coherent state {-d — 0), and 
(ii) S = 4.474160 and x*^' = 0.106O, 



^(3) 

s 
i} 






1.2l7r 
0.7872 



0.5O 


0.457r 
0.9674 



0.4O 
0.45O 
0.527r 
0.9418 



0.3O 

1.06O 

0.47r 

0.9231 



0.2O 

2.10 

0.77r 

0.8751 



0.106O 
4.474160 

TT 

0.9883 



O.IO 
4.80 


0.9924 



0.05O 
9.90 
1.237r 

0.9318 



O.OIO 
49.980 

1.497r 
0.9897 



0.005O 
99.990 

I.Stt 
0.9973 



TABLE II: Field fidelity at f = \U for the field initially in pf 
values of 5 and x*''' with the respective relative phase "d. 



\a){a\ (a 



0) and atom initially in p^ — |e)(e| for different 



which generates a state very close to an odd coherent state 
{•& = tt): those combinations are the only ones that enable the 
generation of an almost Schrodinger cat state at each collapse time 
when the initial field state is either a coherent state or a statistical 
mixture of two coherent states. In appendix B we (analytically) 
show how such results could be understood. 



IV. CONCLUSIONS 

In this work we have investigated the dynamics of a field in a 
lossless cavity interacting with a two-level atom in the presence 
of a nonlinear Kerr-like medium. Using the density operator for- 
malism, we have obtained the exact (RWA) evolution operator 
for this model. We have found that the dynamics of the JCM 
with a Kerr-like medium is considerably richer than shown in the 
literature. The parameters d and x may combine in a way that 



new and interesting features are revealed: for instance, we may 
obtain a periodic dynamics, contrarily to what happens in the 
ordinary JCM. 

The field interacting just with a two-level atom evolves to a 
state not that close to a Schrodinger cat. On the other hand, a 
Kerr-like medium alone may split a single coherent state into a 
superposition of two coherent states, but a very large nonlinear- 
ity is required for that. In any case the generation of pure states 
from a mixed state is not verified. In our model both a Kerr- 
like medium and an atom conveniently detuned from the field are 
important in the nonclassical field generation process. The fine 
tuning of S is important if one wants to match a specific value 
of the nonlinear parameter x i in order to achieve a periodic 
dynamics. Particularly, with a finite detuning there is no need of 
large Kerr nonlinearities, given that in a certain range of parame- 
ters, the larger the detuning the smaller the nonlinearity required. 
However, the detuning does not play a critical role in the gener- 



ation scheme, because even for a zero detuning we have a value 
for the nonhnear parameter (large, though) that gives us a peri- 
odic dynamics. Although schemes based on non-resonant (large 
detuning) interactions do not require precise control over the de- 
tuning, they normally rely upon conditioned measurements. In 
particular, we observed that the periodic dynamics, dictated by 
5^, allows us to recover the initial state at each revival time. We 
have also found that an initially coherent field becomes an almost 
exact superposition of coherent states {Schrodinger cat state) at 
each collapse time. The advantage of our method is that gen- 
eration of sharp Schrodinger cat-like states of the field may be 
achieved in a non-conditioned manner, i.e. without the need of 
collapsing the atomic state. It seems that the non-linearity of a 
Kerr-like medium has an effect on the field that closely resembles 
the non-linear behaviour resulting from an atomic measurement, 
in the case of a conditional method. 

We have also found the conditions in which a field initially in 
a statistical mixture of two coherent states evolves towards an 
almost pure state, e.g. close to an even coherent state. Such 
a quantum field 'purification' may be well understood from the 
phase space point of view: as the field undergoes periodic evo- 
lution, at certain times the overlap of Q-functions coming from 
opposite branches is almost perfect, meaning that an almost pure 
state {Schrodinger cat) has been generated. In order to better 
illustrate the field state generation process, we have calculated 
the cavity field Q-function for successive (close enough) times in 
a way that we could produce an animation of the complete evo- 
lution of the Q-function from i = until half of the revival time. 
The generation of an almost pure field state from a mixed state 
may of course be achieved at the expense of the purity of the 
atomic states, i.e. the atom itself ends up in a mixed state. 
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APPENDIX A: ANALYTICAL RESULTS FOR THE 
DISPERSIVE LIMIT 

We consider here the usual procedure to obtain the dispersive 
Hamiltonian for the JCM 39], but including a Kerr- like medium. 
The dressed states maintain their usual form 

+,n) — sin6'„+i|e,n) + cosO„^i\g,n + 1), (Ala) 

| — ,n) = — cos6„_|_i|e,n) + sinO„^i\g,n + 1), (Alb) 

but the coefficients are given by 



for any 'relevant' n . Under that condition, equation IIA3l l be- 



smy„+i = 



cos 9„ 



\Ji^n + l " In + lf + ^n + 1 

(n„+i — 7„+l) 



(A2a) 
(A2b) 



where uj„+i — 2Q\/n + 1. The corresponding eigenvalues are 

S±,„ = hiooin + i) + hx^'^^n^ ± ifin„+i. (A3) 

The dispersive limit is obtained when we consider i/int, equa- 
tion l|2b), as a small perturbation of the whole Hamiltonian 39]. 
It is equivalent to make 



hQ^ 



i5±,„ ^ hcuoin + i) ± ^n\5\ + nx''^^n{n T 1) ± -r^{n + 1),(A5) 
meaning that we can employ the following effective Hamiltonian 

(A6) 







Analogously to the calculation of section II, we have the evolution 
operator for the dispersive limit given by 



t/f(t)=e--<^>("-")*h"^^'""^" J 
\ e^T' 

so that we may write the evolved field state as 



\i>iii)) = ^c„e 



-iX<^*"^t ix'^'nt -i-r"t 



n), 



(A7) 



(A8) 



where c„ = exp(— |a|^/2)a"/vn! is the coherent state coefficient 
in the number state basis. We are now able to demonstrate the 
initial field state being recovered at each revival time and the 
generation of a Schrodinger cat state at each collapse time: 

Case 1: When 5 = 5^, e.g. 6 = 49.980 and x^^^ = 0.010, we 
recover the initial state at tr = i^/x 



We remark that this result can be easily generalized for any initial 
field state. 

Case 2: Similarly, at ifr — 7r/2x the field evolves to 



IV'f(K))=e-il"l'5]4T( 



\n). 



(AlO) 



If we use e * 2 " 
j), we obtain 



i(l + z)(e- 



) , after multiplying by -75 ( 1 - 



IVf(itr)) = -^(l-ie '5.5xCTQ)-j|je '^^a)) 
v2 

= 7|(l^") -«i-2")). 

in agreement to the numerical result e"* — e*^*^'^ « — i. 



(All) 



APPENDIX B: FROM A STATISTICAL MIXTURE 
TO THE SCHRODINGER CAT 

In the numerical analysis we have noted that the field initially 
in a coherent state with cf) = Q evolves to a Schrodinger cat with 
i9 and the field initially in a coherent state with <p = tt evolves to 
a Schrodinger cat with — 19. Therefore, it is reasonable to suppose 
that the field initially prepared in a statistical mixture of those 
two coherent states evolves to the state 



P( 



i(|5;^){5;tf| + |5;-^){5;-^|). 



(Bl) 



\S\ > UJ„+1, 



(A4) 



-f By relevant we consider the states with significant probability P„ 
(njpjln) of population for the field under consideration. 



i.e. a statistical mixture of two Schrodinger cat states with the 
same relative phase except for a minus sign. Finally, the reason 
why only even and odd coherent states are obtained during the 
evolution of an initial statistical mixture of two coherent states 
becomes clear by noting that 

pf = j5; i9)(5; i?| = C[|5)(5| + |-5)(-5| 
+ cosi9(|5)(-5| + |-5)(5|)-isim?(|5)(-5| - |-5)(5j)], (B2) 



and 



pP-"* = |(|5;^)(5;^| + |5;-^)(5;-^l) 

= C[|5)(5| + |-5)(-5| + cos'i?(|5)(-5| + |-5)(5|)],(B3) 

are equal if and only if i9 = fcyr (fc integer), i.e. only for the even 
and odd coherent states. 
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